Truncated states obtained by iteration 
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Quantum states of the electromagnetic field are of considerable importance, finding potential 
application in various areas of physics, as diverse as solid state physics, quantum communication 
and cosmology. In this paper we introduce the concept of truncated states obtained via iterative 
processes (TSI) and study its statistical features, making an analogy with dynamical systems theory 
(DST). As a specific example, we have studied TSI for the doubling and the logistic functions, which 
are standard functions in studying chaos. TSI for both the doubling and logistic functions exhibit 
certain similar patterns when their statistical features are compared from the point of view of DST. 
A general method to engineer TSI in the running- wave domain is employed, which includes the 
errors due to the nonidealities of detectors and photocounts. 

PACS numbers: 42.50.Dv, 42.50.-p 



I. INTRODUCTION 

Quantum state engineering is an area of growing im- 
portance in quantum optics, its relevance lying mainly in 
the potential applications in other areas of physics, such 
as quantum teleportation [m , quantum computation 2 1 , 
quantum communication Q, quantum cryptography 4|, 
quantum lithography Q , decoherence of states , and 
so on. To give a few examples of their usefulness and 
relevance, quantum states arise in the study of quan- 
tum decoherence effects in mesoscopic fields 0; entan- 
gled states and quantum correlations [1]; interference in 
phase space ; collapses and revivals of atomic inversion 
[lO| |: engineering of (quantum state) reservoirs [T]| : etc. 
Also, it is worth mentioning the importance of the statis- 
tical properties of one state in determining some relevant 
properties of another as well as the use of specific 
quantum states as input to engineer a desired state p^ . 

Dynamical Systems Theory (DST), on the other hand, 
is a completely different area of study, whose interest lies 
mainly in nonlinear phenomena, the source of chaotic 
phenomena. DST groups several approaches to the study 
of chaos, involving Lyapunov exponent, fractal dimen- 
sion, bifurcation, and symbolic dynamics among other 
elements [l^ . Recently, other approaches have been con- 
sidered, such as information dynamics and entropic chaos 
degree |15|. 

The purpose of this paper is twofold: to introduce 
novel states of electromagnetic fields, namely truncated 
states having coefficients obtained via iterative process 
(TSI), and to study chaos phenomena using standard 
techniques from quantum optics, making an analogy with 
DST. We note that, unlike previous states studied in the 
literature [lB|, each coefficient of the TSI is obtained 
from the previous one by iteration of a function. Fea- 
tures of this state are studied by analyzing several of its 
statistical properties in different regimes (chaotic versus 



nonchaotic) according to DST, and, for some iterating 
functions, we found properties of TSI very sensitive (re- 
sembling chaos) to the first coefficient Cq, which is used 
as a seed to obtain the remaining C„. 

This paper is organized as follows. In section 2 we in- 
troduce the TSI and in section 3 we analyze the behavior 
of some of its properties as the Hilbert space dimension is 
increased. In section 4 we show how to engineer the TSI 
in the running-wave field domain, and the corresponding 
engineering fidelity is studied in section 5. In section 6 
we present our conclusions. 



II. TRUNCATED STATES OBTAINED VIA 
ITERATION (TSI) 



We define TSI as 
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n=0 



'Corresponding author: |wesIeybcardoso@gmail.com| 



where C„ is the normalized complex coefficient obtained 
as the nth iteration of a previously given generating 
function. For example, given Co,C„ can be the nth it- 
erate of the quadratic functions: Cnip) = C*^-! + A*; 
sine functions: Cn{^i) — Aisin(C„_i); logistic functions 
CnifJ') = MC'n-i(l — Cn-i); exponential functions: C„ — 
/iexp(Cn-i); doubling function defined on the interval 
[0,1): Cn — 2C„_i mod 1, and so on, /i being a param- 
eter. It is worth recalling that all the functions in the 
above list are familiar to researchers in the field of dy- 
namical systems theory (DST). For example, for some 
values of /i, it is known that some of these functions can 
behave in quite a chaotic manner Also, note that 

by computing all the C„ we are in fact determining the 
orbit of a given function, and because the C„ and P„ , the 
photon number distribution, are related by P„ = |C„p, 
fixed or periodic points of a function will correspond to 
fixed or periodic Pn- Rather than studying all the func- 
tions listed in this section, we will focus on the doubling 
function and the logistic function. These two functions 



2 



have been widely used to understand chaos in nature. 
As we shall see in the following, although very different 
from each other, these functions give rise to different TSI 
having similar patterns. 

III. STATISTICAL PROPERTIES OF TSI USING 
THE DOUBLING AND THE LOGISTIC 
FUNCTIONS 

A. Photon Number Distribution 

Since the expansion of TSI is known in the number 
state |n), we have 

Pn - |C„P. (2) 

Figs. 1 and 2 show the plots of the photon-number dis- 
tribution Pn versus n for TSI using the doubling func- 
tion. The Hilbert space dimension is TV = 50. In order 
to illuminate the behavior of TSI for different values of 
Co, we take Co as 0.3 and 0.29711, respectively shown in 
Figs. 1 and 2. Note the regular behavior for Co — 0.3 
and rather an irregular, or chaotic, behavior for Co — 
0.29711. Figs. 3 and 4 show P„ for the logistic func- 
tion. For Co = 0.2 and /i — 3.49 the logistic function 
behaves regularly (Fig. 3), showing clearly (as in the case 
of the doubling- function) four values for P„; by contrast, 
for /i = 4 and Co — 0.2, Pn oscillates quite irregularly 
(Fig. 4). This is so because the photon number distri- 
bution is equivalent to the orbit of the TSI dynamics 
[l^ . Thus, once a fixed - attracting or periodic - point 
is attained, the subsequent coefficients, and hence the 
subsequent Pn, will behave in a regular manner. Con- 
versely, when no fixed point exists, P„ will oscillate in 
a chaotic manner. Therefore, by choosing suitable Co 
and/or fi, we can compare the properties of TSI when 
different regimes (chaotic versus nonchaotic) in the DST 
sense are encountered. Note the similarity between the 
properties of the logistic and the doubling functions when 
the DST regimes are the same. Interestingly, these simi- 
larities are observed when other properties are analyzed, 
as we shall see in the following. 

B. Even and Odd Photon number distribution 

The functions Podd and Peven represent the photon 
number distribution for n odd and even, respectively, 
give n by Eq. It is well established in quantum optics 
|17| that if Podd > 0.5 the Glauber-Sudarshan P- function 
assumes negative values, prohibited in the usual proba- 
bility distribution function, and the quantum state has 
no classical analog. Since Podd + Peven — 1, the same is 
true when Peven < 0.5 . Figs. 5 and 6 show the behavior 
of Podd for Co = 0.3 and Co = 0.29711 for the doubhng 
function the Hilbert space N is increased. Figs. 7 and 8 
refer to the the logistic function for /i — 3.49 and /i = 4. 
In Figs. 5 and 7 (corresponding to a nonchaotic regime 



0.08- 
0.06- 




Hilbert space dimension N 



FIG. 1: Photon number distribution for the doubling function 
with Co = 0.3. Note a four-period type for the probabilities; 
this regular behavior coincides with nonchaotic behavior of 
the doubling function in the DST sense. 
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FIG. 2: Photon number distribution for the doubling function 
with Co ~ 0.29711. This irregular behavior for P„ coincides 
with the chaotic behavior of the doubling function in the DST 
sense. 



in DST), note that TSI has a classical analog as N in- 
creases. From Figs. 6 and 8 (corresponding to a chaotic 
regime in DST), TSI can behave as a nonclassical state, 
depending on N. More interestingly, note the following 
pattern: whenever the coefficients of TSI correspond to 
the nonchaotic regime in DST, Podd (and so Peven) will 
remain above or below 0.5 on a nearly monotonic curve, 
as seen in Figs. 5 and 7; whenever the coefficients of TSI 
correspond to the chaotic regime in DST, Podd (and so 
Peven) wiU tend to oscillate around 0.5 (Figs. 6 and 8). 



3 




■20 30 40 



50 



FIG. 3: Photon number distribution for the logistic function 
with Co = 0.2 and /i — 3.49. This regular or "four-period" 
type behavior for the probabilities coincides with nonchaotic 
behavior in the DST sense. 
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FIG. 5: Even (solid) and odd (dots) photon number distri- 
butions for the doubling function. We have used Co = 0.3 to 
coincide with nonchaotic behavior in the DST sense. 
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FIG. 6: Even (solid) and odd (dots) photon number distribu- 
tions for the doubling function. Co ~ 0.29711 was chosen to 
coincide with chaotic behavior in the DST sense. 



FIG. 4: Photon number distribution for the logistic function 
with Co = 0.2 and /x = 4. This irregular behavior of P„ 
coincides with the chaotic behavior of the logistic function in 
the DST sense. 

C. Average number and variance 

The average number (n) and the variance (An) in TSI 
are obtained straightforwardly from 



Fig. 9 shows the plot of (n) and Fig. 10 the plot of 
(An) as functions of the dimension N of Hilbert space, 
for the doubling function. Note the near linear behav- 
ior of the average photon number and its variance as 
N increases for Co = 0.3 (nonchaotic regime in DST); 
this is not seen when Co = 0.29711 (chaotic regime in 
DST). Figs. 11 and 12 for the logistic function show es- 
sentially the same behavior when these two DST regimes 
are shown together. 
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and 



(«) = E^W"' (3) 

n=0 D. Mandel parameter and second order correlation 

function 



(An) = ^(n2) - {hf. 



The Mandel Q parameter is defined as 
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FIG. 7: Even (solid) and odd (dots) photon number distribu- 
tions for the logistic function. Here we chose Co = 0.2 and 
= 3.49 to coincide with nonchaotic behavior in the DST 
sense. 



FIG. 9: Average photon number for the doubling function. 
Here we chose Co ~ 0.3 (dots) and Co ~ 0.29711 (solid), 
allowing the Hilbert space A'' to increase to 50. 
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FIG. 8: Even (solid) and odd (dots) photon number distribu- 
tions for the logistic function. Here we chose Co = 0.2 and 
fi = 4 to coincide with chaotic behavior in the DST sense. 



FIG. 10: Variance of the photon number for the doubling 
function. Here we chose Co = 0.3 (dots) as well as Co = 
0.29711(solid), allowing the Hilbert space A'' to increase to 
50. 



(An2 _ (n)) 



while the second order correlation function ^'^•'(O) is 



and for Q < (Q > 0) the state is said to be sub- 
Poissonian (super-Poissonian). Also, the Q parameter 
and the second order correlation function g^^^ are related 
by [3 



If g*^^-* (0) < 0, then the Glauber-Sudarshan P-function 
assumes negative values, outside the range of the usual 

(4) probability distribution function. Moreover, by Eq. ([5|) 
it is readily seen that ^^^■'(O) < 1 implies Q < 0. As 
for a coherent state Q = 0, a given state is said to be a 
"classical" one if Q > 0. 

Figs. 13 to 16 show the plots of the Q parameter and 
the correlation function g^^-'(O) versus N , for both the 
doubling and the logistic functions. Note that TSI is 
predominantly super-Poissonian for these two functions 
{Q > and > 1), thus being a "classical" state 

in this sense for N > 12, while for small values of N 
{N < 12), the Q parameter is less than 0, showing 
sub-Poissonian statistics and is thus associated with a 
"quantum state". From Figs. 13 and 15, note that us- 

(5) ing Co = 0.3 for the doubling function, Cq = 0.2 and 
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FIG. 11: Average photon number for the logistic function. 
Here we chose Co = 0.2 with /li = 3.49 (dots) or = 4 (solid), 
allowing the Hilbert space N to increase to 50. 



FIG. 13: Q parameter for the doubling function. Here we 
chose Co = 0.3 (dots) and Co = 0.29711 (solid), allowing 
the Hilbert space N to increase to 50. Note the irregular 

behavior of the Q parameter when Co = 0.29711, coinciding 
with chaotic behavior in the DST sense. 
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FIG. 12: Variance of the photon number for the logistic func- 
tion. Here we chose Co = 0.2 with /u = 3.49 (dots) or = 4 
(solid), allowing the Hilbert space N to increase to 50. 
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FIG. 14: Second order correlation function for the doubling 
function. Here we chose Co = 0.3 (dots) and Co = 0.29711 
(solid), allowing the Hilbert space N to increase to 50. Note 
the irregular rise of the (;(0) function when Co = 0.29711, 
coinciding with chaotic behavior in the DST sense. 



\x = 3.49 for the logistic function (nonchaotic regime in 

DST), Q shows a linear dependence on N. However, us- 
ing Co = 0.29711 for the doubling function, Cq — 0.2 and 
/i = 4 for the logistic function (chaotic regime in DST), 
Q oscillates irregularly. Similarly, from Figs. 14 and 16, 
using Co = 0.3 for the doubling function, Co = 0.2 and 
fj, = 3.49 for the logistic function, we see that ^''^^(O) 
increases smoothly, while using Co = 0.29711 for the 
doubling function, Co = 0.2 and /x = 4 for the logistic 
function, the rise of g^^^ (0) is rather irregular. This pat- 
tern is observed for other Co as input as well as for other 
values of the parameter jj,, and whenever the dynamics is 
chaotic (regular), the Q parameter and the .g^^-'(O) corre- 
lation function oscillate irregularly (regularly). 



E. Quadrature and variance 

Quadrature operators are defined as 
Xi = |(a + at) ; 

X2=^(a-at). (6) 

where a (a^) is the annihilation (creation) operator in 
Fock space. Quantum effects arise when the variance of 
one of the two quadratures attains a value AXi < 0.5, 
i = 1,2. Figs. 17 and 18 show the plots of quadrature 
variance AXi versus N. Note in this figures that vari- 
ances increase when N is increased. 
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FIG. 15: Q parameter for the logistic function. Here we chose 
Co = 0.2, with = 3.49 (dots) and jj. — 4 (solid), allowing the 
Hilbert space A*' to increase up to 50. Note the irregular be- 
havior of the Q parameter when p = 4, exactly when chaotic 
behavior occurs according to DST. 



FIG. 17: Averaged quadratures for the doubling function 
when the Hilbert space A'' is increased to 50. Doted (dashed- 
doted) line refer to quadrature variance 1 (2) for Co = 0.3; 
solid (dashed) line refer to quadrature variance 1 (2) for 
Co = 0.29711. 
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FIG. 16: Second order correlation function for the logistic 
function. Here we chose Co = 0.2, with fi — 3.49 (dots) and 
fj, = 4 (solid), allowing the Hilbert space A*' to increase up to 
50. Note the irregular rise of g^{0) when ^ = 4, exactly when 
chaotic behavior occurs according to DST. 



FIG. 18: Averaged quadratures for the logistic function with 
Co ~ 0.2 when the Hilbert space N is increased to 50. Doted 
(dashed-doted) line refer to quadrature variance 1 (2) for /i — 
3.49; solid (dashed) refer to quadrature variance 1 (2) for 
M = 4. 



F. Husimi -Q function 



The Husimi Q-function for TSI is given by 



i\TSI 



-mTsi)\\ 



(7) 



where /3 is a coherent state. Figs. 19 to 22 show the 
Husimi Q-function for both the doubhng and the logistic 
functions for N — 15. For the doubling function, we use 
Co = 0.29711 and Co = 0.3, respectively, and for the 
logistic function we use /x = 3.49 and = 4, respectively, 
for the chaotic and nonchaotic regimes. Interestingly, 
even when the chaotic and nonchaotic regimes of the DST 



are compared, Husimi Q-functions show essentially no 
difference from each other. 



IV. GENERATION OF TSI 

TSI can be generated in various contexts, as for exam- 
ple trapped ions [IB], cavity QED "iS, 20], and travelling 



wave-fields |21[ . But due to severe limitation imposed by 
coherence loss and damping, we will employ the scheme 
introduced by Dakna et al. [2l| in the realm of run- 
ning wave field. For brevity, the present application only 
shows the relevant steps of Ref. 21], where the reader 
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FIG. 19: Husimi Q function for doubling function. Here Co 
0.3. 



FIG. 21: Husimi Q function for logistic function. Here Co 
0.2 and = 3.49. 





FIG. 20: Husimi Q function for doubling function. Here Co ~ 
0.29711. 



will find more details. In this scheme, a desired state j^*) 
composed of a finite number of Fock states \n) can be 
written as 



JV 



c 



N 



N 



n(«^-/3;)|o) 



N 



N 



II D{(3k)a+D{f3k) |0) , 



(8) 



k=l 



FIG. 22: Husimi Q function for logistic function. Here Co = 
0.2 and ^ = 4. 



where D{(3n) stands for the displacement operator and 
the Pn are the roots of the polynomial equation 



N 



J2 = 0. 



(9) 



n=0 



According to the experimental setup shown in the Fig.l 
of Ref. 21], we have (assuming 0-photon registered in all 
detectors) that the outcome state is 



JV 



1*) ^l[D{ak+i)a+T^D{ak)\0) 



(10) 



fe=r 
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TABLE I: The roots Pi = \l3k\e~''^^k of the characteristic 
polynomial and the displacement parameters = |Qfc|e~*'''°fc 
are given for TSI using the doubling function for Co = 0.3 
(coinciding with nonchaotic behavior in the DST sense) , = 
5 and T = 0.862. The probability of producing the state is 
0.22%. 



N 






\ak\ 




1 


2.169 


2.638 


1.187 


-0.220 


2 


2.169 


-2.638 


1.155 


1.570 


3 


0.545 


3.141 


1.096 


-2.483 


4 


1.460 


1.084 


1.323 


-2.331 


5 


1.460 


-1.084 


2.225 


1.570 


6 






1.460 


-1.084 



TABLE IIL The roots /^^ = |/3fe|e"**'ffc of the characteristic 
polynomial and the displacement parameters al = |afc|e~*'^°fc 
are given for TSI using the logistic function for Co = 0.2, 
= 3.49 (coinciding with nonchaotic behavior in the DST 
sense), A'^ = 5 and T = 0.893. The probability of producing 
the state is 0.11% 



N 






|Qfc| 




1 


3.948 


3.141 


2.794 


0.051 


2 


0.609 


2.566 


2.195 


-3.045 


3 


0.609 


-2.566 


0.472 


1.570 


4 


1.828 


1.373 


1.830 


-1.959 


5 


1.828 


-1.373 


3.202 


1.570 


6 






1.828 


-1.373 



TABLE II: The roots = |/3fc |e~*'^'3<= of the characteristic 
polynomial and the displacement parameters = lofele"*"^"*: 
are given for a TSI using the doubling function for Co = 
0.29711 (coinciding with chaotic behavior in the DST sense), 
N = 5 and T = 0.867. The probability of producing the state 
is 0.21% 



N 






M 




1 


2.306 


2.692 


1.372 


-0.198 


2 


2.306 


-2.692 


1.130 


1.570 


3 


0.543 


3.141 


1.193 


-2.563 


4 


1.489 


1.089 


1.357 


-2.321 


5 


1.489 


-1.089 


2.289 


1.570 


6 






1.489 


-1.089 



TABLE IV: The roots (31 = \l3k\e~''^i^k of the characteristic 
polynomial and the displacement parameters al = |afc le""^"*: 
are given for TSI using the logistic function for Co ~ 0.2, 
fi = 4 (coinciding with chaotic behavior in the DST sense), 
N — 5 and T — 0.879. The probability of producing the state 
is 0.15% 



N 


\Pk\ 




M 




1 


3.290 


3.141 


2.027 


0.094 


2 


0.563 


2.708 


1.665 


-3.056 


3 


0.563 


-2.708 


0.321 


1.570 


4 


1.893 


1.255 


1.787 


-2.064 


5 


1.893 


-1.255 


3.165 


1.570 


6 






1.893 


-1.255 



where T is the transmittancc of the beam splitter and 
ttfc are experimental parameters. After some algebra, 
the Eq.® and Eq. fTU)) can be connected. In this 
way, one shows that they become identical when ai — 
~T.i=iT-'ai+i and au = T*^-''+\(3k-i - /3fc) for 
k — 2,3, A ... N. In the present case the coefficients C„ 
are given by those of the TSI. The roots = |/3fe| e-"^"" 
of the characteristic polynomial in Eq.(l9]) and the dis- 
placement parameters = \ak \ e~^'^''k are shown in the 
Tables I;II;III and IV, for = 5. 

For N = 5, the best probability of producing TSI is 
0.22% when the doubhng function is used, and 0.15% 
when the logistic funcion is used. The beam-splitter 
transmittancc which optimizes this probability is around 
T = 0.878. 



now take into account the quantum efficiency rj at the 
photodetectors. For this purpose, we use the Langevin 
operator technique as introduced in [22] to obtain the 
fidelity to get the TSI. 

Output operators accounting for the detection of a 
given field a reaching the detectors are given by (2^ 

Sout = ^/?7"m + La, (11) 

where rj stands for the efficiency of the detector and La , 
acting on the environment states, is the noise or Langevin 
operator associated with losses into the detectors placed 
in the path of modes a = a,b. We assume that the detec- 
tors couple neither different modes a, b nor the Langevin 
operators La , so the following commutation relations are 
readily obtained from Eq. lfTTj) : 



V. FIDELITY OF GENERATION OF TSI 

Until now we have assumed all detectors and beam- 
splitters as ideal. Although very good beam-splitters are 
available by advanced technology, the same is not true 
for photo-detectors in the optical domain. Thus, let us 



La, Li 



La, L 



= 0. 



(12) 
(13) 



The ground-state expectation values for pairs of Langevin 
operators are 



9 



? ft 



T ft 



0, 



(14) 
(15) 



which are useful relations specially for optical frequen- 
cies, when the state of the environment can be very well 
approximated by the vacuum state, even for room tem- 
perature. 

Let us now apply the scheme of the Ref.[21] to the 
present case. For simplicity we will assume all detectors 
having high efhciency {rj > 0.9). This assumption al- 
lows us to simplify the resulting expression by neglecting 
terms of order higher than (1 — 77)^. When we do that, 
instead of \TSI), we find the (mixed) state |4'fe) de- 
scribing the field plus environment, the latter being due 
to losses coming from the nonunit efficiency detectors. 
We have, 

X D{aN~i)...a^r'D{ai)\0)Ll 
+ R^-^D{aN+i)a^T^D{aN)a)T^ 
X D{aN^^)...L\r'D{a^)\(S) 

+ i?^-ii?(a^+i)atT"i?(ajv)2]v_i7^" 

X D{aN^i) . ■ ■a)T^D{ai)\Q) 

+ R^'^D{aN+i)%r'D{aN)a^T'' 

X D{aN_i)...a''T^D{ai)\0), (16) 



where, for brevity, we have omitted the kets correspond- 
ing to the environment. Here R is the reflectance of the 
beam splitter, lJ = 1 is the identity operator and Lfc, 
k = l,2..N stands for losses in the first, second . . . N~th 
detector. Although the i^/s commute with any system 
operator, we have maintained the order above to keep 
clear the set of possibilities for photo absorption: the 
first term, which includes = 1, indicates the proba- 
bility for nonabsorption; the second term, which include 
L^i'', indicates the probability for absorption in the first 
detector; and so on. Note that in case of absorption at 
the "k-th detector, the annihilation operator a is replaced 
by the L\, creation Langevin operator. Other possibili- 
ties such as absorption in more than one detector lead to 
a probability of order lesser than (1 — ry)^, which will be 
neglected. 

Next, we have to compute the fidelity [23|, F = 
||(^ , where \^) is the ideal state given by 

Eg. pH)) . here corresponding to the TSI characterized by 
the parameters shown in Tables TIV, and \^fe) is the 
state given in the Eq. fTB]) . Assuming 77 — 0.99, 0.95 
and 0.90 and starting with Co = 0.3 and Co 0.29711 



for the doubling function, we find F ~ 0.9983, 0.9943 
and 0.9909, respectively, and for the logistic function, 
starting with Co — 0.2, ji — 3.49 and = 4, we find 
F ~ 0.9986, 0.9944 and 0.9911, respectively These high 
fidelities show that efhciencies around 0.9 lead to states 
whose degradation due to losses is not so dramatic for 
iV ^ 5. 



VI. COMMENTS AND CONCLUSION 

In this paper we have introduced new states of the 
quantized electromagnetic field, named truncated states 
with probability amplitudes obtained through iteration of 
a function (TSI) . Although TSI can be building using var- 
ious functions such as logistic, sine, exponential functions 
and so on, we have focused our attention on the doubling 
and the logistic functions, which, as is well known from 
dynamical systems theory, can exhibit a chaotic behavior 
in the interval (0,1]. To characterize the TSI for the dou- 
bling and logistic functions we have studied various of its 
features, including some statistical properties, as well as 
the behavior of these features when the dimension N of 
Hilbert space is increased. Interesting, we found a tran- 
sition from sub-poissonian statistics to super-poissonian 
statistics when iV is relatively small {N ^ 12). Besides, 
photon number distribution, which is analogous to con- 
cept of orbits in the study of the dynamic of maps, shows 
a regular or rather a "chaotic" behavior depending on ex- 
isting or not fixed or periodic points in the function to be 
iterated. Interestingly enough, we have found a pattern 
when the properties of TSI for logistic function are com- 
pared with that of TSI for the doubling function from 
the point of view of dynamical systems theory (DST). 
For example, as P„ has an analog with orbits from DST, 
it is straightforward to identify repetitions (or periods) in 
P„, if there are any, when the Hilbert space is increased. 
Surprisingly, although the doubling and the logistic func- 
tion are different from each other, when other properties 
such as even and odd photon number distribution, the 
average number and its variance, the Mandel parameter 
and the second order correlation function were studied, 
they presented the same following pattern: if, from the 
point of view of DST, the coefhcients of TSI for the dou- 
bling and the logistic function correspond to a nonchaotic 
(chaotic) regime, all those properties increases smoothly 
(irregularly) when the Hilbert space is increased. 
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